Magneto-Dielectric phenomena in charge and spin frustrated system of layered iron oxide 
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Dielectric and magnetic phenomena in spin and charge frustrated system M^CU (R is a rare-earth metal 
ion) are studied. An electronic model for charge, spin and orbital degrees in a pair of triangular-lattice planes 
is derived. We analyze this model by utilizing the mean-field approximation and the Monte-Carlo simulation 
in a finite size cluster. A three fold-type charge ordered structure with charge imbalance between the planes 
is stabilized in finite temperatures. This polar charge order is reinforced by spin ordering of Fe ions. This 
novel magneto-dielectric phenomenon is caused by spin frustration and charge-spin coupling in the exchange 
interaction. We show cross-correlation effect in magnetic- and electric-field responses. Oxygen deficiency effect 
as an impurity effect in a frustrated charge-spin coupled system is also examined. 

PACS numbers: 75.80.+q, 72.80.Ga, 75.10.-b, 77.80.-e 



I. INTRODUCTION 

Simultaneous existence of electric and magnetic polariza- 
tions and their control by external field have been recently 
revived as multiferroic phenomena in correlated electron ox- 
ides, ^i 3 ' 4 Behind large coupling between electric and mag- 
netic moments, spin frustration plays dominant roles on mul- 
tiferroic properties. Non-colinear spin structures, such as cy- 
croid, spiral and so on, are realized on a frustrated geometry, 
and spontaneous electric polarization is induced to gain the 
symmetric/anti-symmetric exchange interactions A 6 ^ In this 
viewpoints, this class of materials are recognized as a spin 
driven ferroelectricity. Another class of ferroelectricity is pos- 
sible in correlated electron systems. Apart from the integer 
filling of electron in valence bands, charge degree of freedom 
is active. In particular, around the quarter filling, the long- 
range charge order due to electron correlation is ubiquitously 
observed in several transition-metal compoundsi 8 i 9 i 1Q When 
electronic charge is ordered without inversion symmetry, a 
macroscopic electric polarization appears. This is a ferroelec- 
tricity driven by electronic charge degree of freedom. This 
class of ferroelectricity is realized in low-dimensional organic 
salts, such as the neutral-ionic transition system' 1 '' 2 and a- 
(BEDT-TTF)2l3^ Charge polarized state observed in a lay- 
ered structure manganite Pr(Sro.iCao.9)2Mn207 is attributed 
to the charge-orbital order associated with lattice distortion. 14 
A possibility of ferroelectricty in manganites is also proposed 
in theoretical viewpoint^ 

Rare-earth iron oxides with layered crystal structure 
/?Fe2C»4 (R=Lu, Y, Yb, Er)i& of the present interest belong to 
this class of ferroelectricity. Crystal structure of /?Fe2C<4 con- 
sists of paired Fe-O triangular-lattice layers and R-O blocks 
stacked along the c axis. Schematic view of a paired Fe-O 
layer, termed the W-layer, is shown in Fig. 0I a )- Average 
valence of Fe ions is +2.5, implying that equal amounts of 
Fe 2+ and Fe 3+ occupy the W-layer. Therefore, this mate- 
rial is recognized as a spin-charge frustrated system. Charge 
structure was investigated by the electron and x-ray diffrac- 



tion experiments 
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In LuFe2C>4, below 500K, streak- 



type diffuse scattering was observed along (1 /3 1 /3 /) lines, 
and below 320K, spots appear at (1 /3 1 /3 3m + 1 /2) in the 
streak lines associated with zigzag modulations. In this pa- 



per, we use the hexagonal index, although the space group 
is R3m. These experimental results are interpreted as two- 
and three-dimensional charge orders of electrons. The three- 
dimensional order of Fe 2+ and Fe 3+ was also confirmed 
by the resonant x-ray scattering technique at Fe .ST-edge^ 
As for the magnetic properties, magnetization in LuFe2C<4 
starts to increase around 250Ki 2 -L Neutron diffraction exper- 
iments revealed that magnetic Bragg peaks at (1/3 1/3 m) 
appear and a ferrimagnetic order realizes below this temper- 
aturei 22 ' 23 ' 24 ' 25 ' 26 Electric polarization and dielectric anoma- 
lies were observed around the three-dimensional charge or- 
dering temperature, although the dielectric constant shows 
strong dispersive and diffusive naturei 20 ' 27 Several magneto- 
dielectric phenomena were also reported around the ferri- 
magnetic ordering temperaturei 20 i 28 ' 29 It is worth noting that 
these dielectric and magnetic phenomena depend on the rare- 
earth metal element, R, and the oxygen stoichiometry; in 
YFe204, with decreasing temperature, the three fold-type 
charge order is changed into a four fold-type one which is 
extremely sensitive to oxygen deficiency i 30 ' 3 'i 32 These micro- 
scopic and macroscopic experiments denote that 3d electronic 
charges are responsible for the dielectric anomalies, and cou- 
ple strongly with spins. 

To elucidate mechanism of dielectric phenomena in 
/?Fe204, Yamada and coworkers proposed a model for the 
three fold-type charge orderi 8 This charge-structure model is 
shown in Fig.[5ja) which will be introduced in more detail in 
Sect. [TITl This is a \/3 x \/3 structure in a plane, and along 
[110], electronic charges are aligned ■ ■ -Fe 3+ Fe 3+ Fe 2+ - ■ ■ in 
the lower plane and ■ • Fe 3+ Fe 2+ Fe 2+ - ■ ■ in the upper one. 
That is, electronic charge is polarized between the upper and 
lower planes, and finite electric dipole moments exist in the 
W-layer. Based on this polar charge model and the neutron 
diffraction data, a possible spin structure in the ferrimagnetic 
ordered phase was proposed^ 

A number of the experimental result o 16,18 ' 20 ' 21 ' 22 and the 
analy se a 17 i 24 - 33 i 34 i 35 imply that electronic processes and in- 
teractions are crucial for novel dielectric properties in this 
material. In this paper, we present a microscopic theory 
of electronic structure and magneto-dielectric phenomena in 
/?Fe2C>4. We focus on 3d electronic structure in a W-layer 
which is a minimum and main stage for the low-energy elec- 
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FIG. 1: (a) : A pair of triangular-lattice planes (W-layer). (b) : A 
FeC>5 cluster. 



tronic state. We first suggest the orbital degree of freedom 
in a Fe 2+ ion, and derive an electronic Hamiltonian in a W- 
layer. This model consists of the long-range Coulomb inter- 
actions and the exchange interaction derived from the gener- 
alized pd model. We analyze the charge structure by using the 
mean-field approximation and the Monte-Carlo (MC) simula- 
tion. The three-fold type polar charge order competes with 
other type non-polar ones, and is stabilized at finite tempera- 
ture. This is caused by charge fluctuation in a frustrated tri- 
angular lattice. We furthermore examine spin structure and 
coupling between spin ordering and electric polarization. The 
polar charge order is strongly stabilized below the magnetic 
ordering temperature. This magneto-dielectric phenomenon is 
attributed to spin frustration in a triangular lattice. We demon- 
strate novel electric and magnetic responses which are avail- 
able to examine the present theoretical scenario. Effects of 
oxygen deficiency on electric polarization are also studied. 

In Sect, mi we derive the model Hamiltonian for electronic 
structure in a W-layer. In Sect. HID numerical results for 
the charge structure and electric polarization are presented. 
Calculated results for the spin structure and the magneto- 
dielectric responses are shown in Sect. [TV] Examined oxy- 
gen deficiency effects are introduced in Sect. V. Section[VT]is 
devoted to discussion and concluding remarks. Preliminar 
results for the present study have been published in Refs 
andHH Study of a doubly degenerate orbital model in a hon- 
eycomb lattice as an orbital model for RF^Oa is presented in 
separate papers! 34 ' 36 



generate orbitals (—ad^ + bd^^, ad yz + bd xy ) with E', and 
(ad x 2_ y 2 + bd ZXl — ad xy + bd yz ) with E". Numerical coeffi- 
cients a and b satisfy the relation a 2 + b 2 = 1 . We obtain 
that the degenerate E' orbitals take the lowest energy with 
b = 0.89, and the first excited level is E". The energy dif- 
ference between E' and E", AE e i_ e h, is about O.leV which 
is smaller than that between E" and A', AE E n_ A i ~ 0.6eV. 
When we see the crystal structure in detail, an Fe ion is not 
located at center of a O5 cage. Distance between the Fe ion 
and the O3 plane denoted by h [see Fig.[TJ(b)] is about 0.1 A 
in LuFe2C>4 i 37 ' 38 We obtain that, with taking h into account, 
AE E /_ E « increases and AE E n_ A i decreases. The hybridization 
effects between Fe 3d and O 2p orbitals may increase these 
level separations. However, because of the small value of 
A£ e /_ e h, we do not exclude a possibility that (d x 2_ y 2, d xy ) 
and (d u , d yz ) couple strongly with each other, i.e. a ~ b, and 
that the E" level is the lowest. In any cases, the lowest orbitals 
are degenerate. As a result, in Fe 3+ , each orbital is singly 
occupied, and total spin 5 = 5/2 of the high-spin state. On 
the other hand, in Fe 2+ , one of the degenerate lowest levels is 
doubly occupied by a hole, and 5 = 2. Thus, two fold orbital 
degeneracy exists in Fe 2+ . This is represented by the orbital 
pseudo-spin operator defined by 



(1) 



HH's 



where dj^ s is the creation operator for an Fe 3d hole with or- 
bital Ll, spin s(=t, I) at site i, and <7^/ is the Pauli matrices. 
In following part of this paper, we assume for simplicity that 
the two orbitals in the lowest level are (d x 2_ y 2,d xy ), and the 
index 11 in Eq. (TJ) takes the two. The z component of the op- 
erator Tf is 1/2 (—1/2) for the state where a hole occupies 
the d x 2_ y 2 (dxy) orbital. Even in the case where the orbitals 
in the lowest level are (d zx ,d yz ), the following part of this pa- 
per is valid by reinterpreting that the index ji in Eq. (TJi takes 
{dyz, d zx ). 

We set up the model Hamiltonian for the electronic struc- 
ture in a W-layer. The 3d electrons in the Fe ions and the 2p 
ones in O which hybridizes with Fe 3d are introduced. We 
start from the following generalized pd Hamiltonian, 



II. MODEL HAMILTONIAN 

We start from the electronic structure in a single Fe ion in 
the W-layer. This ion is five-fold coordinate with three O ions 
in the xy plane and two at apices as shown in Fig. [JJ(b). We 
calculate the crystalline-field splitting of the Fe 3d orbitals in 
the FeOs cluster. Five O ions are replaced by point charges 
with valence of — 2e, and their positions are determined by 
the crystal structure data i 37 ' 38 The hydrogen-like wave func- 
tions are adopted for the Fe 3d orbitals, and the effective 
nuclear charge is taken to be +8. The split 3d orbitals are 
identified by the irreducible representation in the group: 
the dj )Z 2_ r 2 orbital with A', and two sets of the doubly de- 



■^pd = + + ' + ffi/ , (2) 

with 

-\ L ^dJuAtic'd^^din'a, (3) 



3 




FIG. 3: Inter-site Coulomb interactions between Fe ions. Solid, 
broken and dotted arrows represent interactions between the near- 
FIG. 2: Three two-dimensional coordinates, (a x , ce y ), (j3 A ,/3 v ) and est neighbor (V cNN ), the next nearest neighbor (V a bNN) and the third 
(Yx, Y y ) in a FeO triangular lattice. Filled and dotted circles represent neighbor (V cNNN ) Fe-Fe bonds, respectively. 
Fe and O ions, respectively. 
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ions are taken into account in the last term of Eq. (f2]), i.e. 
j£y . We consider the largest three interactions in the W-layer 
as shown in Fig. [3] the inter-plane NN interaction (V c nnX 
the intra-plane NN one (V a bNN) an d the inter-plane next NN 
one (VcNNn)- This is because (1) these Coulomb interac- 
tions are a minimum set which reproduces the three-types of 
charge structures observed experimentally, and (2) a distance 
between the 4th neighbor Fe ions in the W-layer is compara- 
ble to that between the NN W-layers. This will be discussed 
in Sect. Hn] in more detail. Summations in J%v take the three 
kinds of pairs. When the 1 /r-type Coulomb interaction is as- 
sumed, we obtain V c nn /VabNN = 1.2 and V c nnn /VabNN = 0.77 
for LuFe2C<4. By introducing the pseudo-spin operator Qj for 
charge degree of freedom, J%y is rewritten as an antiferromag- 
netic Ising model 



where rfj is the creation operator for the Fe 3d hole with or- 
bital p. (= xy, x 2 — y 2 , yz, zx, 3z 2 — r 2 ) and spin a (=1,1) 
at site ;', and p\ yg is for the O 2p hole with orbital V (= 
x, y, z). Number operators are defined by nf^ = dj^di^a, 

n 'jvo = P'jvoPjv°> and n i = Lno^fiia- A simbo1 S n is a con " 
necting vector between Fe and NN O ions along directin rj. 

Interactions in a Fe ion are described in the 1 st term of Eq. ©, 
Jtfj, where the level energy e^, the intra-orbital Coulomb in- 
teraction U d , the inter-orbital one W d , and the exchange in- 
teraction I d are considered. Interactions in J4? p are defined in 
the same way with those in J^/. Hopping of a hole between 
the nearest neighboring (NN) Fe and O ions in the same plane 
is described in Jtf with the transfer integral t pd . For conve- 
nience, we introduce the three two-dimensional coordinates 
(rf x ,riy) with 77 = (o£,j3,y), which are obtained by a rotation 
of the xy axis by 2nm^j3 with (m a: mp : my) = (0,1,2) (see 
Fig.|2]i. In each coordinate, we define the operators as 
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(8) 



In the bond direction 77, the c/ 7) 2_ J) 2 and p^. orbitals compose 
the rj bond. The inter-site Coulomb interactions between Fe 



(cNN abNN cNNN \ 

E E VabNN + E VcNNN Q$&j, (9) 

W (y) M / 

where a constant term is omitted. The operator Qf takes 1/2 
and —1/2 for Fe 3+ and Fe 2+ , respectively. The charge con- 
servation is imposed by a relation £ ; Q : - = 0. 

Based on the extended pd Hamiltonian J%? pc i, we derive the 
effective Hamiltonian for the superexchange interactions be- 
tween NN Fe ions in a plane. This interaction arises from 
virtual hopping of holes between Fe ions. The Hamiltonian 
is derived by the 4th order projection-perturbation procedure 
in terms of the hopping term Jf pc i- Following two exchange 
processes are considered: 



d M p Q d N ->d M+l p°d N - 1 



d M p°d N , 



and 



d M p°d N 



d M - l p 2 d N - l ^d M p°d N , 



(10) 



(11) 



where we adopt the hole picture, and M and N represent 
the numbers of holes. These are termed the dd- and dpd- 
processes, respectively, from now on. Here, we present the 
outline of derivation, and details are given in Appendix A. A 
general form of the Hamiltonian is 
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Q&tP, (12) 



where P is the projection operator for the d 5 (d 4 ) high-spin 
states in Fe 3+ (Fe 2+ ), Q = 1 — P, and £, is the initial-state en- 
ergy. Many body effects of My in the intermediate states are 
considered approximately; we assume that the intermediate- 
state energies for d M+l p°d N ~ l [see Eq.dTObl are higher than 
the initial- and final-state energies by a constant energy pa- 
rameter V, which is of the order of the inter-site Coulomb in- 
teraction. Then, we set Mq = Md + M p + V. We interpret that 
effects of J^y in other intermediate states are included in the 
charge-transfer energy, Act, between the d ri 2_ rj 2 and p^ y or- 
bitals. The obtained Hamiltonian Mj is classified by valences 
of Fe ions in the initial states, i.e. Fe m+ -Fe" + where n and to 
take 2 and 3, and electron configurations in the intermediate 
states denoted by k. The Hamiltonian is given as 



•%j= E 



)-* 




(13) 



FIG. 4: The lowest-energy spin and orbital configurations for a 
Fe 2+ -Fe 2+ bond (a), for Fe 3+ -Fe 2+ (b), and for Fe 3+ -Fe 3+ (c). 
Open, filled and dotted circles represent Fe 2+ , Fe 3+ and O ions, re- 
spectively. Spin and orbital configuration in (a) is energetically close 
to ferromagnetic spin alignment with d x 2_ y 2 and dp a orbitals. 



All possible intermediate states k are taken into account in 
Eq. (13) which consists of 6 terms in M^ 11 ^ and M {23) - k , 
and 4 terms in M^~ k . Explicit formulae of all terms are 
presented in Appendix A. Here, we show some representative 
terms: 



^(22)- 



7( 22 )- 1 £(I ; .I 7 + 6)fi-2T ; - 7) ,T^.) 

«/> v y 



-G5 



(14) 



Jf(23)-1 = ,/(23)-l 



E (Jrlj + 

(ij) V 



2 Xjnj 



(15) 



^(33)-l =7 (33)-l^/ J .. J ._^ 

(ij) V 4 ' 



(16) 



We define the spin operators I, and J, for Fe 2+ and Fe 3+ with 
amplitudes of 2 and 5/2, respectively. The orbital operator is 
redefined in the (r\ x , rj y ) coordinate as 



Zir, = Tf cos ( ym, J+Tj sin I 



(17) 



This operator takes 1/2 (—1/2), when the d rj 2_ ri 2 {dr\ x ^ y ) 
orbital is occupied by a hole. In a given pair of i and j 
sites, subscripts 77, and rjj in t (J)i . and T/^. are automati- 
cally determined. The exchange constants are defined by 

y( 23)-l = _ 2 ,2 dc/[25A(23) _ i]; 

and y' 33 ' -1 = 4/J dc /[25A(33\_ 1 ] where f^/ c is the transfer inte- 

(^COS0)/ACT 

We introduce 



gral between NN Fe ions defined by tddc 
with the Fe-O-Fe bond angle 6(= 120° 



the intermediate-state energies as A 



*(23)-l 



V, and A 



(33)-l 



(22)- 1 " 

4/'' + y. 



W d - l d + V, 
It is worth to 



note that 1) JfK 22 )~' is expressed as a product of charge, spin 
and orbital interactions between given sites i and j, and 2) 
j^(32)-/ i nc i uc i es a linear term of the orbital pseudo spin be- 
cause Fe 3+ dose not have the orbital degree of freedom. 

After all, we obtain the Coulomb- and exchange-interaction 
Hamiltonian 



(18) 



where 3%y and M 3 ] are given in Eqs. (0 and ( fT3l , respectively. 
Before going to the numerical results calculated in the Hamil- 
tonian, we briefly mention the energy scales of charge, spin 
and orbital degrees of freedom, and signs of the exchange 
interactions. The inter-site Coulomb interactions provide a 
larger energy scale than the exchange interactions. Thus, the 
charge sector is frozen at the highest temperature in compar- 
ison with spin and orbital ones. This is consistent with the 
experimental results in LuFe2C<4 where the charge ordering 
temperature (about 320K) is higher than the spin ordering one 
(about 250K)] 17 ' 18 i 24 By calculating the exchange energy in a 
given NN bond, we estimate stable spin and orbital config- 
urations. This is not trivial from the Goodenough-Kanamori 
rule because of the 120° bond angle. The energy parameter 
sets for the exchange coupling constant /('""W are determined 
from the experimental data in LaFe03 , 39 i 40 and I, and J; are 
assumed to be Ising spins. We obtain the spin and orbital con- 
figurations for the lowest exchange energies as 1) for a Fe 2+ - 
Fe 2+ bond, Iffj = —4 (antiferromagnetic) and X\ = Ty = 1 /2, 
which is energetically close to Iffj = 4 (ferromagnetic) and 
T; = -Tj = 1 j2, 2) for Fe 3+ -Fe 3+ , JfJj = -25/4 (antiferro- 
magnetic), and 3) forFe 2+ -Fe 3+ , IfJj — 5 (ferromagnetic) and 
Tj = 1/2. Schematic views for the stable configurations are 
presented in Fig.|U In the neutron scattering experiments, the 
ferrimagnetic phase indexed as (1/3 1/3 to) appears around 
250K. Possible magnetic structures are shown in Fig. [13] 
which will be explained in more detail later. In this structure, 
Fe 2+ ions in the 2Fe 3+ -Fe 2+ (upper) plane are surrounded by 
six NN Fe 3+ . Thus, the exchange Hamiltonian in this plane 
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is reduced into a form of L(; ; -)[(l/2) ± T^.] which becomes 
a constant by using the relation of Y,(if) T ir/; = 0. This rela- 
tion is also applicable to the Fe 2+ -Fe 3+ bonds in the 2Fe 2+ - 
Fe 3+ (lower) plane where three Fe 2+ -Fe 3+ bonds connecting 
a Fe 2+ ion are equivalent [see Fig. [T3l . Therefore, the orbital 
part of the exchange Hamiltonian in this ferrimagnetic phase 
is mapped onto the following orbital model defined on a Fe 2+ 
sublattice; 



orb 



(19) 



(a) OOifi [010] (b) co 1/4 




••cct* 

<X>«KKJ 



where (c a ,ep,ey) represent the three unit vectors connecting 
NN Fe 2+ sites in a honeycomb lattice. A summation Y!i takes 
Fe 2+ sites in one of the two sublattices in a honeycomb lattice. 
The coupling constant J 0I \, is given by the exchange constants 
j{nm)-k^ j n tn j s mo( jel, j{ j s obtained theoretically that the 
orbital does not show a conventional long-range order down 
to very low temperature of the order of 0.005/ or b. Therefore, 
for simplicity, we assume that the pseudo-spin operators for 
orbital in ,3tfj are set to be zero in the following calculation. 
Theoretical study of the orbital model on a honeycomb lattice 
is presented in separate papers! 34 ' 36 
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III. CHARGE STRUCTURE AND ELECTRIC 
POLARIZATION 

First we focus on the charge structure and the electric po- 
larization by analyzing the inter-site Coulomb interaction term 
J%v. We apply, at the first stage, the mean-field approximation 
to J%v, and obtain stable charge structures. The charge con- 
servation is taken into account by adding the chemical poten- 
tial term, — VextLiGJ, i n the Hamiltonian. We assume that the 
expectation value (<2?) is periodic along the (110) or (210) di- 
rections, and takes the same amplitude along (110) or (010), 
respectively. Periodicity L is taken up to 12. In upper and 
lower planes, (2f)'s are independent and have the same pe- 
riodicity along the (110) or (210) directions. Each solution 
is characterized by the momentum q = (M /2L,M /2L,n) or 
(M/2L,0,«) where M is the number of nodes of (Q--) along 
the (110) or (210) directions, respectively. When a phase 
difference between (<2f)'s in the upper and lower planes is 
(7r), n takes (1 /2). Phase diagram is determined by com- 
paring the free energy. Representative charge structures are 
shown in Fig. [5] Four types of CO's in this figure, denoted 
by CO 1/3, CO 1/4, CO 1/2 -I and CO1/2-II, are characterized 
by momenta q = (1/3, 1 /3,0) = q 1/3 , (1 /4, 1/4, 1 /2) = q 1/4 , 
(1/2,1/2,0) = q 1/2 _i and (1/2,0,0) = q 1/2 _ n , respectively. 
As suggested by Yamada and coworkers, CO 1/3 shows finite 
electric polarization due to charge imbalance between the two 
triangular-lattice planes i 17 ' 18 A ratio of Fe 2+ and Fe 3+ is 1 : 2 
(2 : 1) in the upper (lower) plane. In other charge structures, 
equal numbers of Fe 2+ and Fe 3+ occupy the upper and lower 
planes, and there is no electric polarization. 

Mean-field phase diagram at zero temperature is presented 
in Fig. [6] The non-polar CO^-H and CO 1/4 structures are 
stable in the regions of Vcnn /Vcnnn < 2 and Vcnn /Vcnnn > 



FIG. 5: Charge structures in a W-layer : (a) CO 1/3, (b) CO 1/4, (c) 
CO1/2-I, and (d) CO1/2-II. Filled and open circles represent Fe 3+ 
and Fe 2+ , respectively, and large and small circles are for Fe ions in 
the upper and lower planes, respectively. Lower panel in each figure 
is a side view from [110]. 



2, respectively. The polar CO 1/3 structure appears only on the 
phase boundary where CO 1/3 is degenerate with COi /2-II and 
CO1/4. Realistic parameter values for RF^Oa correspond to 
a shaded area in Fig. [6] We fix a value of V C NN /VabNN to be 
1.2, as shown by a dashed line in Fig. [(J and calculate finite- 
temperature phase diagram (Fig. [7]). The polar COj n is stabi- 
lized in a wide region between CO 1/4 and the CO 1/2 -II. 

Beyond the mean-field calculation, we examine the charge 
structure in finite temperature by using the MC simulation. 
To avoid a trap of a simulation in local minima, we adopt the 
multi-canonical MC (MUMC) method, 41 Simulations are per- 
formed on a paired triangular lattice of L x L x 2 (= 2N) (L = 6 
and 12) sites with the periodic -boundary condition in the ab 
plane. We use 6 x 10 6 MC steps to obtain a histogram in the 
MUMC method and 16 x 10 6 MC steps for measurement. We 
calculate the charge correlation function and the electric po- 
larization P defined by 



N(q) 
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(2A0 : 
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with 



p={p 2 yi\ 



(20) 



(21) 



(22) 
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FIG. 6: Mean-field phase diagram for charge order at zero tem- 
perature. Three charge structures, CO1/4, CO1/2-II, an d CO1/3, are 
degenerate on a line of V c nn = 2K c nnn, and the two, COj/4 and 
CO; m-I, are degenerate on a line of V c nnn = 0. Shaded area corre- 
sponds to a region for i?Fe2C>4. Phase diagram in finite temperatures, 
shown in Fig. [7] is calculated on broken line. 
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FIG. 7: Mean-field phase diagram for charge structure in finite tem- 
peratures. The Coulomb interaction VcNn/Vann is chosen to be 1.2. 



where r,- is a position of site ;, and represents a summa- 
tion of site i in the upper (lower) plane. 

The charge correlation functions at V c nnn /VabNN = 0.58, 
0.6 and 0.62 are presented in Fig. [8] At V c nnn /V a bNN = 
0.6, Af(q!/ 3 ) shows a hump around T = r/VabNN = 0.18 
and keeps a finite value down to the low temperature limit. 
From the specific heat data, we identify T = 0.18 corresponds 
to the charge-ordering temperature. In the both cases of 
VcNNN/VabNN = 0.58 and 0.62, A^qm) is dominant in high 
temperatures, starts to decrease around T = 0.2 and disap- 
pears at the lowest temperature. On the contrary, the charge 
correlation A^(q 1 /2-n) and A^q^) grow up around T = 0.05, 
and increase with decreasing temperature. These results and 
the specific heat data imply that the charge order at qi/3 is 
changed into the other type of charge order at qi/4 (qi/2-1) 
around f =0.04 (0.045) for V cNNN /V abN N = 0.58 (0.62). Tem- 



FIG. 8: Charge correlation functions N(q) at V C NNN/V a bNN = 0.60 
(a), 0.58 (b) and 0.62 (c) calculated in ,ffiy. The Coulomb interaction 
VcNN /VabNN is chosen to be 1.2. 
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1 ' v abNN 

FIG. 9: Electric polarization Pas a functions of V c nnn/ VabNN cal- 
culated in My . The Coulomb interaction V c nn /VabNN is chosen to 
be 1.2. 



perature dependence of P at several values of V c nnn /VabNN is 
presented in Fig. [9] At V C NNN/V a bNN = 0.6, P remains down 
to the low temperature limit. Apart from V c nnn /V a bNN = 0.6, 
P starts to decrease at the temperature where Af(qi/2-n) and 
iV(qi/4) grow up, and disappears at the lowest temperature. 
These results obtained by the MUMC method are qualitatively 
consistent with the ones in the mean-field calculation. 

The polar charge structure characterized by qi/3 and the 
transition to the another structure characterized by q^u at 
Vcnnn /VabNN < 0.6 are consistent with the experimental re- 
sults. In LuFe2C>4, charge order indexed as (1/3 1/3 3m + 
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FIG. 10: Two-sublattice structure in the CO1/3 phase. Amplitude of 
the mean-field at Fe sites surrounded by broken circles is 0.9V a bNN 
and that at other sites is 2. 1 V abNN . 

1/2) appears around 350K and remains, at least, down to 
around 20K£L On the other hand, in YFe2C>4, charge order 
indexed as (1/3 1/3 3m + 1 /2) observed at room temperature 
is changed into the one as (1/4 1/4 3/4) around 250K^iH 
We suppose that different rare-earth metal ions slightly change 
ratio of the Coulomb potentials, and LuFe2C>4 (YFe204) 
corresponds to the parameter region of V c nnn /V a bNN ~ 0-6 
(VcNNN/VabNN < 0.6) in the present calculation. 

Stability of the CO 1/3 phase is attributed to large thermal 
fluctuation in the geometrically frustrated lattice. A key is- 
sue is the two-sublattice structure in this charge ordered phase 
(see Fig. [TOlii^ 3 - Fe 2+ ions in the 2Fe 2+ -Fe 3+ (lower) plane 
and Fe 3+ in the Fe 2+ -2Fe 3+ (upper) one belong to a sublattice 
termed A. Other Fe ions belong to another sublattice termed 
B. All in-plane NN ions of a site in the sublattice B have an 
opposite valence. On the other hand, a site on the sublattice 
A is surrounded by three NN Fe 2+ and three NN Fe 3+ in the 
plane. Therefore, the Coulomb potentials at these sites from 
the in-plane NN ions are canceled out, and charge fluctuation 
is able to occur easily without loss of V a bNN- It is obtained 
in the numerical calculation that an amplitude of the mean- 
field on the sublattice A is 0.9V a bNN at low temperature which 
is much less than that on the sublattice B, 2.1V a bNN- Large 
charge fluctuation at the sites grows up with increasing tem- 
perature, and contributes to the entropy gain at finite tempera- 
ture. On the contrary, in the CO ^2 -II and CO ^4 structures, all 
Fe 2+ (Fe 3+ ) are equivalent and charge fluctuation is weaker 
than that in the sublattice A of CO 1/3. This is the reason why 
the polar charge order characterized by (1/3, 1/3, 0) is more 
stable than other charge structures in finite temperatures. 

Let us focus on the charge structure in low temperatures 
in more detail. As shown in Figs.[8](b) and [9] saturated val- 
ues of A^qi/3) andP in V C NNN/VabNN=0.6 at the low tempera- 
ture limit are 0.032 and 0.094, respectively, which are smaller 
than the values expected from the ideal CO^ phase, 0.056 
and 0.33 respectively. This implies that the charge config- 
uration at low temperature in Vcnnn/Vann = 0.6 is not the 
ideal CO1/3 state. We analyze the probability histogram in 
the MUMC simulation, and examine the charge configura- 
tions realized in lowest temperatures. These are classified 
into the following three configurations: the polar CO1/3 struc- 
ture shown in Fig. [5] (a), partially polarized charge structures 



characterized by the momentum qi/3, termed COa, and non- 
polar structures termed COb- Detailed structures of COa an d 
COb are shown in Appendix IE1 In COa, the polarization is 
P =N /3—n\/~N with an integer number n satisfying a relation 
of < n < /3. Degeneracy of a sum of these configura- 
tions is of the order of 2^/3^" ~ 2^. As for the COb 
state, degeneracy of the configuration is also of the order of 
Because of the coexistence of these charge structures, 
the saturated values of P and A^q^) are smaller than the ex- 
pected values from the ideal CO1/3 state. This tendency is 
remarkable in the large system size, as shown in Fig. [9] 

This coexistence of the polar and non-polar CO states im- 
plies that the full polarization expected from the ideal CO 1/3 
state is realized by additional weak interactions. The long- 
range Coulomb interactions between the NN W-layers is 
one of the candidates. This scenario is plausible, since, in 
LuFe204, the electric polarization appears around the three- 
dimensional charge-ordering temperature ! 20,27 We examine 
effects of the inter W-layer Coulomb interaction based on a 
model where two W-layers stacked along the c axis are cou- 
pled by the Coulomb interaction. Saturated values of N(<\\p) 
and P at low temperature are identical to the expected val- 
ues from the ideal polar CO1/3 state. Roles of the exchange 
interaction as another candidate to lift the degeneracy are ex- 
amined in the next section. 



IV. SPIN STRUCTURE AND MAGNETO-ELECTRIC 
EFFECT 

In this section, we introduce spin degree of freedom and 
examine coupling between the electric polarization and the 
magnetic ordering. The Hamiltonian j£y + Mj is analyzed 
by utilizing the MUMC method in a 6x6x2-site cluster. 
The spin operators I, and J, in J^j are assumed to be Ising 
spins because of the strong magnetic anisotropy observed in 
/?Fe204. 24 The energy parameters in the Hamiltonian are cho- 
sen to be U d = 7.8, W d = 6.2, I d = 0.8, U p = 4.1, W p = 
2.9, P = 0.6, tP d = 1.8, Act = 3 and V = 1 in a unit of 
VabNN- These are determined from the experimental date in 
LaFe03 i 39 i 40 In this section, the orbital pseudo-spin operators 
in Jtj are set to be zero, as explained in Sect. [II] In particular, 
we focus on a parameter region around V C NNN/VabNN = 0.6, 
where CO1/3 is seen down to the lowest temperature in Fig. [7] 
and that around 0.58-0.59, where the transition from CO1/3 to 
COi u is shown in Fig. [7] 

Temperature dependences of the charge correlation func- 
tion, the spin correlation function defined by 

%) = 7^2l^) e "'' q ' (ri " rj0 . ( 23 ) 

where Kf = If (Jf) for Fe 2+ (Fe 3+ ), and the electric polariza- 
tion are calculated. Results at V c nnn/ VabNN = 0.60 and 0.59 
are shown in Figs. fTTI and fT2l respectively. For comparison, 
we also plot the data obtained in My. At V c nnn /VabNN = 0.6, 
three characteristic temperatures, T — 0.2, 0.085 and 0.015, 
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FIG. 11: Charge correlation functions N(q) (a), spin correlation 
functions S(q) (b), and electric polarization P (c) calculated in + 
Mj. Dashed lines in (a) and (c) are results obtained in J£y. Parame- 
ters are chosen to be V C NN/VabNN = 1-2 and V c nnn /^abNN = 0.60. 



are seen in A^q^). The highest one, 0.2[= 7co(qi/3)], cor- 
responds to the charge ordering temperature for CO 1/3. Other 
two are the magnetic ordering ones at momentum qi/3. At 

f = 0.085 [= r N (q 1/3 )] and 0.015, spins in the Fe 2 +-2Fe 3+ 
and 2Fe 2+ -Fe 3+ planes in CO1/3 (see Fig. [3]l start to order, 
respectively. This double-magnetic transition may be an ar- 
tifact in the present model where the inter-plane exchange 
interactions are neglected, and spins in the upper and lower 
planes are independent with each other. We expect that the 
inter-plane exchange interactions are much smaller than the 
in-plane ones. This is because, when electrons in the dxy and 
d x 2_ y 2 orbitals are concerned, there are no exchange paths 
between Fe ions in an inter-plane NN bond. When higher- 
order exchange processes and/or contributions from other d 
orbitals are taken into account, weak inter-plane interactions 
may unify the double transition in the present calculation. As 
shown in Fig.QT| the charge correlation function at qi/3 and 

the polarization increase at T = 0.085 and 0.014. Results 
clearly show that magnetic ordering enhances stability of the 
polar COj/3 phase. In the low temperature limit, ^(qi/3) and 
P take 0.056 and 0.33 respectively, which are the ideal values 
in CO^. At VcNNN/VabNN = 0.59 (Fig. [T2j, a weak shoul- 
der in iV(qi/3) around T = 0.2 corresponds to the charge or- 
dering for COj/3. Sequential charge ordering transition oc- 
curs from CO1/3 to CO1/4 around T = 0.015 [= 7co(qi/4)], 




T/V abNN 



FIG. 12: Charge correlation functions Af(q) (a), spin correlation 
functions 5(q) (b), and electric polarization P (c) calculated in Jtfy + 
Jtfj. Dashed lines in (a) and (c) are results obtained in Parame- 
ters are chosen to be V C NN/VabNN = 1-2 and V C NNN/V a bNN = 0.59. 



which is lower a little than the result in JZy. Magnetic order at 
qi/3 appears around f = 0.1 [= 7/v(q 1/3 )]. Below 7co(qi/4). 
magnetic structure is also changed; the spin correlation func- 
tions atqi/4 and (5/12, 5/12, 0) become dominant. It is also 
shown, in this parameter, that the electric polarization is en- 
hanced in the CO 1/3 phase. A similar temperature dependence 
is obtained in V a bNN /^cNNN = 0.61, where the CO1/2-H phase 
appears in low temperatures instead of CO 1/4. 

Low temperature charge and spin structures at 
VcNNN/VabNN=0-6 are shown in Fig.QJla). Charge structure 
is identified to be CO1/3. Spins at Fe 3+ in the Fe 2+ -2Fe 3+ 
(upper) plane and those at Fe 2+ in the 2Fe 2+ -Fe 3+ (lower) 
one are aligned antiferromagnetically. On the other hand, 
spin directions at Fe 2+ in the Fe 2+ -2Fe 3+ plane and at Fe 3+ 
in the 2Fe 2+ -Fe 3+ one are not determined uniquely. We note 
that the spin structure in the Fe 3+ -2Fe 2+ (lower) plane is 
sensitive to the parameter values in Jtfj. The structures shown 
in Figs. [iJl a) and (b) are almost degenerate with each other. 
The numerical results presented in this paper are obtained in 
the parameter sets where the spin structure in Fig. \l3[ a) is 
obtained. However, qualitative difference for the results in 
the two parameter sets is not seen. It is also true that essence 
of the coupling between the spin ordering and the electric 
polarization shown in Fig.QT]does not depend on the detailed 
parameter values. Since the antiferromagnetic interaction 
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FIG. 13: Charge and spin structures in the polar COj/3 phase at 

^cNNN AabNN = 0-6. Filled and open circles represent Fe 3+ and 
Fe 2+ , respectively, and large and small circles are for Fe ions in the 
upper and lower planes, respectively. Arrows represent spin direc- 
tions. At Fe sites surrounded by broken circles, spin directions are 
not determined uniquely. Spin structures in (a) and (b) are almost 
degenerate. 



(a) C0 1/4 (6) C0 1/2 - II 




FIG. 14: Charge and spin structure at V^nn/V^nn = *- a '' anc ' 
that at 0.61 (b). 

between NN Fe 3+ -Fe 3+ bonds in the 2Fe 3+ -Fe 2+ (upper) 
plane is robust, Fe 2+ spins are surrounded by three up and 
three down spins in their NN Fe 3+ sites. Therefore, spin 
directions in Fe 2+ are not determined uniquely as explained 
above. Because the number of these sites is N/3, there 
is a macroscopic number of degenerate spin states of the 
order of 2 N ^ which contributes to the entropy gain in finite 
temperatures. This is a kind of partially disordered phase, 
which has been examined in the antiferromagnetic Ising 
model on a triangular latticed In the present case, spins in 
Fe 2+ and Fe 3+ are inequivalent, i.e. 5 = 2 and 5/2, and this 
partial disordered state becomes more stable in comparison 
with that in the conventional Ising model. Since this spin 
structure is realized in the CO 1/3 structure and the spin 
entropy is larger than the charge entropy in the non-polar 
and partially-polar charge ordered phases, i.e. COa and 
COb, the polar CO 1/3 is reinforced through the spin-charge 
coupling in the exchange Hamiltonian. This is a kind of 
"order from fluctuation" mechanism, and, in the present 
spin-charge coupled system, a ferroelectric order is stabilized 
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FIG. 15: Magnetic-field effect in charge correlation function N(q) 
(a), and electric polarization P (b). Parameters are chosen to be 
VcNN/VabNN = 1-2 and VcNNN/VabNN = 0.6. Inset of Jb) shows 
magnetic-field dependence of the electric polarization at T =0.05. 



by spin fluctuation. This phenomenon is not expected in 
COj/2-I, COj/2-II and CO^. Low temperature charge and 
spin structures in V c nnn /VabNN = 0.59 and 0.61 are shown 
in Fig. [14] In both cases, all spins in NN Fe 2+ -Fe 2+ and 
Fe 3+ -Fe 3+ bonds are aligned antiferromagnetically. There 
are a number of degenerate spin states; for example, when 
all spins on a chain along [110] in COm-II are flipped, 
the exchange energy is not changed. However, this spin 
degeneracy is of the order of 2^, which is smaller than 
0(2^) in C0 1/3 . 

In a remaining part of this section, we examine responses 
to the electric and magnetic fields in the present novel spin- 
charge coupled system. First we pay our attention to the 
magnetic-field effect by introducing the Zeeman term of the 
Hamiltonian 

Mk=H^Kf, (24) 

i 

where Kj = If or Jf for Fe 2+ or Fe 3+ , respectively, and H is 
the magnetic field. The Hamiltonian Jty + ,Jtfj + J#h is ana- 
lyzed by utilizing the MUMC method in a 6 x 6 x 2-site cluster. 
Magnetic field dependence of the electric polarization and the 
charge correlation functions at V c nnn /VabNN = 0.6 and 0.59 
are presented in Figs. [15] and [T6J respectively. Temperature 
in Fig. [15] is chosen to be T =0.05, which is below the Neel 
temperature 7n(Qi/3)> ar, d those in Fig. [16] are T =0.05 and 
0.01, which are between 7n(Qi/3) an d tne charge ordering 
temperature of COw 4 [7co(Qi/4)]> an d below 7co(<li/4)> re- 
spectively. When VabNN is taken to be leV, magnetic field 
ff/VabNN = 0.01 corresponds to about lOOTesla. In the mag- 
netically ordered CO1/3 phases at VdsnsrN/VabNN^O.6 ar, d 0.59, 
applying the magnetic field reduces the electric polarization. 
On the other hand, in the antiferromagnetic CO ^4 phase at 
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FIG. 16: Magnetic-field effect in charge correlation function N(q) 
(a), and electric polarization P (b). Parameters are chosen to be 
VcNN/VabNN = I- 2 and V cNNN /V ab NN = 0-59. Inset of (b) shows 
magnetic-field dependence of the electric polarization. Solid circles 
and squares are calculated at T =0.01 and 0.05, respectively. 



FIG. 17: Electric-field effect of spin correlation function 5(q) at 
qi/3. Parameters are chosen to be V c nnn /^abNN = 0.6 (a) and 0.58 
(b) with VcNN/^abNN = 1.2. Insets show electric-field dependence of 
^(ll/3)- Temperatures are T =0.05 in (a), and 0.03 in (b). 



V C NNN/VabNN = 0.59 [see T < 7fco(qi/4) in Fig. [Mb)], the 
electric polarization is induced by applying the magnetic field. 
At the same time, the charge correlation function iV(q 1/3) in- 
creases andAf(qi/4) decreases. Similar results are obtained at 
VcNNn/Kann = 0.61, where the COi/2-n phase collapses and 
the electric polarization appears below T = 0.015 by apply- 
ing the magnetic field. Thus, opposite magnetic-field effects 
are obtained in the magnetically ordered CO 1/3 phase and the 
antiferromagnetic COj 14 and COj n-n- 

We, first, pay our attention to the negative magnetic-field 
effect in the magnetically ordered COi/3 phase. As explained 
in Sect. [HI] the three charge structures, the polar CO1/3, the 
partially polar COa and the non-polar COb, coexist at H = 0. 
Among the three, the polar COj/3 is a dominant structure, be- 
cause of the large spin entropy due to the spin degeneracy of 
the order of 2 N ^. By applying the magnetic field, these N/3 
spins are aligned to be parallel to the magnetic field, and the 
macroscopic spin degeneracy is lifted. On the other hand, in 
both the non-polar COb and the partially polar COa, a macro- 
scopic degeneracy in the charge configuration, which is of or- 
der of 2^, survives under the magnetic field. As the result, 
the charge entropy in COa and COb overcomes the spin one in 
CO1/3, and P is reduced. In other words, under the magnetic 
field, the present spin-charge coupled system is mapped onto 
a spin-less model described by My where the charge entropy 
plays a dominant role. On the contrary, the positive magnetic- 
field effect in the antiferromagnetic CO 1/4 phase is explained 
from the Zeeman energy. Under the magnetic field, the fer- 
rimagnetic structure in the polar CO1/3 phase is more stable 
than the antiferromagnetic one in CO1/4, and the polarization 
appears below 7co(<Ji/4)- However, under a high magnetic 
field larger than ///V a bNN ~ 0.01, the polar CO1/3 competes 



with COa and COb, and the polarization decreases, as dis- 
cussed above. 

We are also able to demonstrate the response to the electric 
field in the present spin-charge coupled system. The static 
electric field E applied along the c axis is formulated by the 
Hamiltonian 

J4? E = -eEd 1 £0, (25) 

i 

where d is a distance between the two triangular-lattice 
planes in a W-layer. Amplitude of the electric field 
eEd/V^nn = 0.01 corresponds to about 50MV/m, when we 
take VabNN =leV and d = 2. 2 A. Electric-field dependence of 
the spin correlation functions at V C NNN/VabNN=0.6 and 0.58 
are presented in Fig. [17] By applying the electric field, 
the spin correlation at qiy 3 is enhanced, in particular, below 
7co(qi/4) in V C NNN/VabNN=0.58. This is a consequence of the 
polar COj/3 phase stabilized by the electric field. The results 
would be use as a test for the present scenario. 

V. EFFECT OF OXYGEN DEFICIENCY 

It is well known that several dielectric and magnetic prop- 
erties in /?Fe204, e.g. charge and spin ordering temperatures, 
are extremely sensitive to the oxygen stoicheometry, denoted 
by 5 in the formula /?Fe204_g . 3 ° l32144 Effects of the oxygen 
deficiency in this system are recognized as the impurity effects 
in charge-spin coupled system with geometrical frustration. 
Here we examine roles of oxygen deficiency on the magneto- 
electric phenomena. We simulate following two aspects of the 
oxygen deficiencies; 1) charge imbalance between Fe 2+ and 
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FIG. 18: Oxygen-deficiency effects of electric polarization P at 
VcNNN/VabNN = 0.6 (a), and 0.58 (b). Bold, broken and dotted lines 
are for result without deficiency, that with charge-imbalance effect, 
and that with electro-static potential, respectively. 

Fe 3+ , which is introduced by the modified charge conserva- 
tion relation as N~ Qf = —25, and 2) random electro-static 
potential around defect sites. This is modeled by the Hamilto- 
nian 

Jft = 2££V R (|i-j|)j25, (26) 

• i 

where and represent summations of defect sites and that 
of the neighboring Fe sites, respectively. We assume that a 
defect site is in the FeO plane, and effective charge of a defect 
is 2+. Amplitudes of the electro-static potentials are estimated 
by the 1 / r-type potential as V R = 1 .73V a bNN an d 1 .60VabNN 
for the NN and NNN sites from a defect site, respectively. 
The model Hamiltonian My + Jffj + Jf%. is analyzed with the 
relation A^ 1 Q] = -25 by the MUMC method. One defect 
is introduced in a 6 x 6 x 2 site cluster. This concentration 
corresponds to 5 = 0.05. 

In FigQjO oxygen-deficiency effect on the electric polariza- 
tion is presented. In CO1/3 [see below 7co(qi/3) in Fig.fTSla) 
and between ?co(qi/3) and Tcofairt) m Fig.fT8tb)l, both the 
charge-imbalance and electro-static potential effects suppress 
the electric polarization. On the contrary, in the antiferromag- 
netic CO 1/4 phase below T^(q l ^ 4 ) in Fig. [T8l b), the electric 
polarization is induced by both the two-types of deficiency ef- 
fects. These results are consistent with the electron-diffraction 
experiments in YF^O^gr^- in samples with large 5, the four 
fold-type charge order disappears, but the three fold-type in- 
dexed as (1/3 1/3 3m + 1 /2) is robust. 

We turn to explain a mechanism of the charge-imbalance 
effect. Reduction of P in the CO1/3 phase shown in Fig.[T8la) 
is a kind of an usual impurity effect which tends to break the 
long-range order. On the contrary, increase of P in low tem- 
peratures shown in Fig. \T%\ b) is related to the characteristic 



charge frustration in CO1/3 as follows. The charge imbalance 
represented by a relation A^ 1 £,<2f = —25 implies replace- 
ment of some Fe 3+ ions in a stoichiometric system by Fe 2+ . 
This corresponds to flipping of pseudo-spins Q^. It is rather 
trivial, in Fig. [5] that this flipping of a Qj happens uniquely 
in the CO1/2-F CO1/2-H and CO1/4 structures. However, in 
COi/ 3 , there are two ways to flip a Q { because of the two- 
sublattice structure mentioned in Sect. [TTTJ Fe 2+ sites sur- 
rounded by NN three Fe 2+ and three Fe 3+ in a plane (sub- 
lattice A), and those surrounded by six NN Fe 3+ in a plane 
(sublattice B). A pseudo-spin in the sublattice A is able to 
be flipped easily in energy. We numerically calculate energy 
costs due to a flipping in sublattice A is about 40% of that in 
sublattice B, and is about 65% in CO1/2-II and CO 1/4. Such 
low-energy charge excitations in CO 1/3 stabilize the charge 
structure under the charge imbalance. 

The electro-static potential effect is also understood from 
a viewpoint of a soft charge structure in the CO1/3 phase. 
Since an effective charge of a defect is 2+, Fe 2+ ions, rather 
than Fe 3+ , tend to assemble to screen this positive excess 
charge. However, due to the Coulombic interaction between 
Fe 2+ ions, a simple cluster consisting of Fe 2+ around a defect 
is not energetically favored. Exchange of Fe 2+ and Fe 3+ be- 
tween the planes in a W-layer is able to reduce such Coulom- 
bic energy. Energy cost for this kind of exchange of Fe 2+ and 
Fe 3+ is much lower in COj n than that in other charge ordered 
structures. That is, the electro-static screening for excesses 
charge easily occurs in CO1/3 because of the two-sublattice 
structure. 



VI. DISCUSSION AND CONCLUDING REMARK 

Here we have remarks on some issues which are not in- 
cluded explicitly in the present model and calculation. Effects 
of the electron transfer in 3d orbitals are not taken into ac- 
count in the Hamiltonian Eq. < fT~8T >. This may be reasonable 
for the first-step theoretical model in /?Fe204. It is because, 
even above the three-dimensional charge-ordering tempera- 
ture (250K) in YFe204, the electric resistivity p shows an 
insulating behavior; p increases with decreasing temperature. 
A magnitude of p about 250K is of the order of 10 2 i2cm, 4 - 
which is much larger that that above the Verwey transition in 
Fe^04^- Therefore, we suppose that dominant electron mo- 
tion is caused by thermal motion, rather than quantum elec- 
tron transfer. This is supported by the experimental data in 
the dielectric constant; it shows strong dispersive feature well 
described by the Debye model based on the thermal fluctu- 
ation of dipole moments i 20 ' 27 ' 47 Electron-transfer effects for 
the charge ordered phase in a triangular lattice have been in- 
vestigated for some low-dimensional organic salts.iMi In the- 
oretical calculations based on the V — t and extended Hubbard 
models at quarter filling, a metallic phase appears in a param- 
eter region between two different charge orders, or it coex- 
ists with the three-fold type charge order. We suppose that 
small electron transfer in RVziOa stabilizes the CO 1/3 phase, 
although diffusive features in the dielectric function becomes 
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more remarkable. 

Lattice degree of freedom and a coupling with electron 
are not included explicitly in the present calculation. In our 
knowledge, there are no detailed crystal structure data in spin- 
charge ordered phases. It is thought from the experimental 
analyses in YFe2C>4 that the crystal symmetries in both the 
two- and three-dimensional charge-ordered phases indexed as 
(1/3 1/3 3m + 1 /2) are trigonal, but that in the four-fold type 
charge order is monoclinic. 44 This result indicates that the lat- 
tice distortion in the three-fold type charge order is weaker 
than that in other phases. This is consistent with the present 
results for a soft charge-order character in CO1/3; amplitude 
of the charge correlation function is smaller than that in other 
phases. A weak lattice distortion expected in the three-fold 
charge order is also related to the symmetry of the CO 1/3 
structure, where the rhombohedral symmetry remains in a 
FeO planes, unlike other charge ordered phases. 

In Sect. IIVI We show that the electric polarization is rein- 
forced by the ferrimagnetic ordering, and this originates from 
the spin entropy due to the frustrated geometry. The long- 
range exchange interactions and/or the magneto-striction ef- 
fects, which are not included explicitly in the present model, 
may release the spin degeneracy. In these cases, we suppose 
that the spins, which are not fixed in the present model [see 
Fig- El, are loosely bounded by such low-energy scale inter- 
actions. However, these still fluctuate in a temperature region 
which is higher than the energy scale of the interactions, and 
contribute to the entropy gain. 

In the present paper, we analyze an electronic model de- 
fined in a single W-layer which is recognized as a minimal 
and main stage in /?Fe204. Obtained results provide a start- 
ing point to elucidate a variety of magnetic and dielectric 
phenomena. We briefly discussed, in Sect. [Ill] some roles 
of the inter W-layer Coulomb interaction. To clarify the 
three-dimensional charge and spin structures ^ 17 ' 25 ! 26 ! 30 and the 
magneto-dielectric response along the c axis, a more realistic 
modeling for the inter W-layer interactions, in particular, the 
inter W-layer exchange interactions, and analyses of a three- 
dimensional model are necessary. 

In summary, electronic structure and magneto-dielectric 
phenomena in the rare-earth iron oxides with frustrated ge- 
ometry are examined. The model Hamiltonian describing the 
electronic interactions between charge, spin and orbital de- 
grees of freedom of Fe ions is derived. This model is analyzed 
by utilizing mainly the Monte-Carlo technique in a finite size 
cluster. The three fold-type charge order associated with elec- 
tric polarization is stabilized in finite temperature in compar- 
ison with two and four fold-type charge structures. The two- 
sublattice structure in this polar charge order plays a crucial 
role. This phase is reinforced by the magnetic ordering, due 
to the spin frustration and the coupling between charge and 
spin in the exchange Hamiltonian. Novel magneto-electric 
responses to the external fields are available as a test of the 
present scenario. Effects of the oxygen deficiency are under- 
stood from the viewpoint of impurity effects in a frustrated 
spin-charge coupled system. Through the present study, we 
provide a unified scenario for a variety of magnetic and di- 
electric phenomena in /?Fe204. 
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FIG. 19: Intermediate states of the exchange processes in a Fe 2+ - 
Fe 2+ bond represented by hole picture, (a), (b) and (c) are for the 
Hamiltonian Jfi 22 )- 1 , ^( 22 )~ 2 , and JT( 22 )~ 3 in the <ta-processes, 
respectively, and (d), (e), and (f) are for ( 22 >~ 4 , ^( 22 ^ 5 , and 
-ones, respectively. Long and short arrows rep- 
resent spins with 5 = 3/2 and 1 /2, respectively. 



APPENDIX A: EXCHANGE HAMILTONIAN 

In this appendix, we show details of the superexchange pro- 
cesses and an explicit form of the Hamiltonian Jffj introduced 
in Sect. HU There are two kinds of the superexchange pro- 
cesses termed the dd- and t/pcZ-processes as introduced in 
Eqs. ( [Tol l and ( TTTb . The Hamiltonian J%j is classified by va- 
lences of Fe ions, Fe" 1+ -Fe" + in the initial and final states, 
and the electronic structure in the intermediate states k [see 
Eq. dl31) l. In this appendix, nearest neighboring Fe sites con- 
cerning in the superexchange interactions are denoted as i and 
j. Electron configuration in Fe and O ions are represented in 
a hole picture. 



1. Exchange Interactions in Fe 2+ -Fe 2+ 

For the aW-processes, electron configurations in the inter- 
mediate exchange processes are denoted as d 3 p°d 5 . Five 
holes are at a site j and three holes at ; with S = 3/2. The 
intermediate states are classified by the spin and orbital states 
at site j [see Fig. [T9l a)-(c)l: (a) the total spin of Fe holes at site 
j, S, is equal to 5/2 and both the two E' orbitals are occupied, 
(b) S= 3/2 and the two E' are occupied, and (c) S = 3/2 and 
one of the E' is occupied. The explicit forms of the exchange 
Hamiltonian are given by 
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Here we define the projection operators for charge 

1 



P Q± = -±O z 
2 U " 



and those for orbital 



pz± _ 



(A4) 



(A5) 



The exchange parameters are given as 

/(22)-l = 

-d/c/[10A (22M ], J {22) - 2 = ^ C /[10A (22) _ 2 ], and 
y (22)-3 = f2 rfc /[4A ( 2 2 )_ 3 ], and A( TOl) _ i is the energy of the sec- 
ond order intermediate states given by A( 22 )_ i = W d —I cl + V, 
A (22) _ 2 = W d + M d + V, and A (22) _ 3 = U d + \I d + V. We 

define t ddc = t] d cos /Act and t dds = t dd sin 0/A C t- 

In the intermediate states of the c/pii-process, two holes oc- 
cupy the O ion. These states are classified by the spin and 
orbital states at the O site, [see Fig. [Hld)-(f)]: (d) the total 
spin of the O holes, S, is equal to 1 and both the p x and p y 
orbitals are occupied by holes, (e) S = and two p orbitals 
are occupied, and (f) S = and one of the p orbitals occupied 
by holes. The exchange Hamiltonians are given by 
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FIG. 20: Intermediate states of the exchange processes in a Fe 2+ - 
Fe 3+ bond represented by hole picture, (a), (b) and (c) are for the 
Hamiltonian Jfi 23 )- 1 , .j^( 23 )- 2 , and JT( 23 )" 3 in the <ta-processes, 
respectively, and (d), (e), and (f) are for ( 23 >~ 4 , ^f 23 ^ 5 , and 
in the dpd-ones, respectively. Long, medium and short 
arrows represent spins with S = 2, 3/2 and 1 /2, respectively. 
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The exchange parameters are given as 

7 (22)-4 = 

-^ s /[4A( 22) - 4 ], J {22) - 5 = d/s/[4A (22 )- 5 ], and 
J {22) ~ 6 = <ifc/[2A ( 22)-6] with A (22) _ 4 = 2A CT + W 
A ( 22 )-5 = 2A CT + W" 7 +IP, and A (22) _ 6 = 2A CT + U p . 



2. Exchange Interactions in Fe 2+ -Fe 3+ 

Electron configurations in the intermediate states are 
d 4 p°d 5 and d 4 p 2 d i for the dd- and afpcZ-processes, respec- 
tively. As well as the exchange interaction in the Fe 2+ -Fe 2+ 
bond, these are classified by the spin and orbital structures in 
the d 5 and p 2 sites for the dd- and lip^-processes, respectively 



(see Fig. 



The explicit forms of the exchange Hamiltoni- 



Jf(32)-1 = /(32)-l 



p^-pQ+pQ- 
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^(32)-6 = y (32)-6 £ (J; . j _ 5)P J+pQ+pQ- j (AM) 

(y) 

for the dpd-ones. The exchange parameters are given as 

' (32M = -2^ c /[25A (32) -i], J^- 2 = 2^ C /[25A (32) _ 2 ], 
^ (32) - 3 = & c /[K)A ( 32)-3], J^ 4 = -^ S /[5A (32) _ 4 ], 
^ (32) ~ 5 = d/ s /[5A( 32 )- 5 ], and/ 32 '- 6 = 2f 2 , c /[5A (32) _ 6 ] with 
5I d + V , A (32) _ 3 = U d - W d + M d + V, 
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(32)-2 : 



A (32) _ 4 = 2A CT + W - A (32) _ 5 = 2A CT + W p + I p , and 



M32)-6 
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B 3+ 



3. Exchange Interactions in Fe + -Fe 

Electron configurations in the intermediate states are 
d 4 p°d 6 and d 4 p 2 d 4 for the dd- and rfpif -processes, respec- 
tively. In the d 6 configuration for the aW-process, total spin is 
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FIG. 22: Schematic pictures of the CO 1/3 structure (a) and one of 
COA(b). When, in CO1/3, Fe 3+ in the upper plane and Fe 2+ in 
lower one on chains indicated by broken lines are exchanged, the 
COa structure in (b) is obtained. 
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FIG. 21: Intermediate states of the exchange processes in a Fe 3+ - 
Fe 3+ bond represented by hole picture, (a) is for the Hamiltonian 
in the aW-processes, respectively, and (b), (c), and (d) are 
for ,if ( 33 >~ 2 , ^( 33 ^ 3 , and ,if( 33 ^ 4 in the dpd-ones, respectively. 
Long, medium and short arrows represent spins with 5 = 2, 3/2 and 
1/2, respectively. 



2 [see Fig. l2Tl a)l. and the explicit form is given by 



jg&y 1 = 7 (33)- 1 £ j j. . j . _ ^ ) pQ+pi 
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FIG. 23: One of the CO B structures. 



The exchange parameter is 7' 33 ) 1 = 4fj (/c /[25A 



i(33)-l. 



with 



A(33)_i = U d +4I d + V. For the c/paf-processes, the interme- 
diate states are classified by the spin and orbital structures in 
the O site [see Fig.|2TJb)-(d)]. The Hamiltonians are given by 



^(33)- 



— I 

4 / ' J 



dQ+ p Q+ 



(A17) 



^(33)-4 =-/ (33)-4£ fj..j._^\pQ+pQ+_ (A18) 

The exchange parameters are /( 33 '~ 2 = — 4f 2 v/s /[25A(33)_ 2 ], 

J (33) - 3 = 4^ s /[25A (33) -3] and = 8t 2 d J [25 A {33) _ 4 ] 

with A (33) _ 2 = 2A CT + W - P, A ( 3 3) _3 = 2A CT + W p + P, 
and A( 3 3)_4 = 2A CT + U p . 



APPENDIX B: CHARGE STRUCTURES OF CO A AND CO B 

In this appendix, detailed charge structures in the COa an d 
COb phases, introduced in Sect. [TTTJ are presented. Charge 
configurations of COa are constructed from CO 1 p . Start from 



the CO1/3 structure shown in Fig. 1221 a). and focus on chains, 
e.g. along [110], where different valences of Fe ions occupy 
the upper and lower planes. Let exchange all Fe 2+ and Fe 3+ 
in any of these chains with each other. One of the obtained 
configurations, termed COa, is shown in Fig. l22l b). These 
structures of COa are energetically degenerate with COip in 
the Hamiltonian J%v + The Coulomb interaction between 
the 2nd NN sites in the plane may lift the degeneracy. When 
the number of the chains, where Fe 2+ and Fe 3+ ions are ex- 
changed, is n (0 < n < 2\//V/3), the electric polarization is 
P = N/3 — n\/N. The degeneracy of a sum of these states is 

of the order of 2^/3^" ~ 2 • Such exchange of Fe ions 
are also allowed on chains along [120] and [210] directions. 

In another degenerated structure, COb, the configuration in 
one side of the W-layer is constructed by stacking two kinds 
of chains alternately. These chains are schematically given as 
• ••OOIOOI- (termed chain and ■■■••o»«o-- - (chain 
S3) along the [110] direction where • and o represent Fe 3+ and 
Fe 2+ , respectively. As shown in Fig. [23] without energy loss 
of VabNN, there are two ways to stack a chain stf on a chain 58, 
and vice versa. These are denoted as A and A, and B and S in 
Fig H3 Therefore, these configurations are degenerated of the 
order of . Charge configuration on another side of the W- 
layer is uniquely determined to gain the inter-plane Coulomb 
interactions V c nn an d Knnn- Obtained charge structures are 
degenerate with the CO 1/3 structure. 
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